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(57)Abstract: 

PROBLEM TO BE SOLVED: To provide a method 
and device for optimizing an illumination profile for a 
pattern of selected patterning means. 
SOLUTION: The method and device for 
microlithography includes optimizing illumination 
modes, based on the features of a specific mask 
pattern. The illumination is optimized by determining 
an appropriate illumination mode, based on the 
diffraction order of a reticle and the autocorrelation of 
a projection optics. By eliminating the part of the 
illumination pattern which has no influence on 
modulation, excess DC light can be reduced, thereby 
improving the depth of focus. Optimization of the 
mask pattern includes addition of sub-resolution 




features, to reduce pitches and truncate the probability density function of a space width. 
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1 Title of Invention 

ILLUMINATION OPTIMIZATION FOR 
SPECIFIC MASK PATTERNS 

2 Claims 

1 . A method of optimizing an illumination profile for a selected patterning 
means pattern, comprising the steps of: 

defining a transmission cross coefficient function for an optical system mcluding 
an illuminator and the selected patterning means pattern; 

determining relative relevance to imaging of diffraction orders based on the 
selected pattern; and 

calculating an optimized illumination configuration from the transmission cross 
coefficient function, weighting regions of the ulumination configuration based on the 
relative relevance to imaging of the diffraction orders. 

2. A method as in claim 1 , wherein the step of determining relative 
relevance to imaging of diffraction orders further comprises the step of: 

deterrnining a characteristic pitch of the selected mask partem. 

3. A method as in claim 2, farther comprising the steps of: 
identifying a critical region of the selected pattern prior to determining the 

characteristic pitch, 

wherein determining the characteristic pitch for the selected partem is performed 
by determining the characteristic pitch for the critical region. 

4. A method as in claim 3, wherein the step of identifying a critical region 
further comprises identifying a plurality of critical regions and wherein determining the 
characteristic pitch for the critical region includes: 

comparing a pitch of each identified critical region; and 

if the pitch of each identified region is substantially equal, determining the 

characteristic pitch for the critical region to be equal to the characteristic pitch of one of 

the identified regions. 



(25) 2002-334836 

5. A method as in any one of the preceding claims, further comprising 
weighting regions of the illuminator configuration based on a selected optimized metric 
selected from the group consisting of depth of focus, end-of-line, image log slope (1LS), 
image slope (IS), and aberration sensitivity. 

6. A method as in claim 1 or 2, further comprising the steps of: 

identifying a plurality of critical regions; 

determining a pitch of each of the identified critical regions; and 

calculating optimized illumination configurations from the transmission cross 

coefficient function, weighing orders based on relevance to imaging of diffraction 

orders for each critical region, and 

wherein calculating an optimized illumination configuration further comprises 

calculating a composite optimized illumination configuration based on the calculated 

optimized illumination configuration for each critical region. 

7. A method as in any one of the preceding claims, further comprising the 
steps of: 

modifying the selected pattern by optical proximity correction techniques to 
reduce a total number of different pitches in the mask pattern. 

8. A method as in claim 7, wherein the step of modifying the selected 
pattern by optical proximity correction techniques further comprises addition of sub- 
resolution features to the selected mask pattern. 

9. A method as in claim 7, wherein the steps of modifying the selected 
pattern and calculating an optimized illumination configuration steps are iterated. 

10. A computer program for optimizing an illumination profile comprising 
program code means that, when executed on a computer system, instruct the computer 
system to cany out the steps of the method of any one of claims 1 to 9. 

11. A device manufacturing method comprising the steps of: 

(a) providing a substrate that is at least partially covered by a layer 
of radiation-sensitive material; 

(b) providing a projection beam of radiation using an ilhimination system; 

(c) using patterning means to endow the projection beam with a pattern 
in its cross-section; 
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(d) projecting the patterned beam of radiation onto a target portion of the layer of 
radiation-sensitive material, wherein, prior to step (d), the cross-sectional intensity 
distribution in the projection beam produced in step (b) is tailored to the pattern 
employed in step (c) using a method according to any one of claims 1 to 9. 

12. A lithographic projection apparatus comprising: 

- an illumination system for providing a projection beam of radiation; 

- a support structure for supporting patterning means, the patterning 
means serving to pattern the projection beam according to a desired 
pattern; 

- a substrate table for holding a substrate; 

- a projection system for projecting the patterned beam onto a target 
portion of the substrate, 

wherein the apparatus additionally comprises; 

- calculating means, for defining a transmission cross coefficient function for the 
ilhiminator and the patterning means, determining relative relevance to imaging of 
diffraction orders based on the pattern produced by the patterning means, and 
calculating an optimized illumination configuration from the transmission cross 
coefficient function, weighing regions of the illumination configuration based on the 
relative relevance to imaging of the diffraction orders; and 

- selecting means, for selecting the cross-sectional intensity 
distribution in the projection beam exiting the illumination system in 
accordance with the illumination configuration calculated by the said 
calculating means. 

1 3. A method of optimizing a selected mask design comprising: 
identifying critical features of the selected mask design; 

determining an optimized illumination profile based on diffraction orders of the 
critical features; and 

modifying the selected mask design by use of optical proximity correction 
techniques which are selected to reduce a number of pitches present in the selected mask 
design. 

14. A method as in claim 1 1 wherein the optical proximity correction Anther 
comprises addition of sub-resolution features selected to modify a continuous 
probability density function of the space width of the selected mask design such that the 
modified probability density function has an increased discretization. 
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15. A method as in claim 13 or 14 wherein said step of determining an 
optimized illumination profile comprises the steps of the method of any one of claims 1 
to 9. 

1 6. A computer program for optimizing a selected mask design comprising 
program code means that, when executed on a computer, instruct the computer to carry 
out the method of any one of claims 13 to 15. 

3 Detailed L^iiipiiar. of Invention 

This invention relates generally to a method and apparatus for microlithographic 
imaging. More particularly, it relates to an apparatus and method for optimizing an 
illumination configuration according to the specific pattern being imaged. 
Optical lithography systems are in current use in the manufacture of integrated circuits 
and odicr fine featured products such as programmable gate arrays. In a most general 
description, a lithography apparatus includes an illumination system which provides a 
projection beam of radiation, a support structure which holds a patterning means, a 
substrate table which holds a substrate, and a projection system (lens) for imaging the 
patterned beam onto a target portion of the substrate. 

The term patterning means should be broadly interpreted as referring to devices 
and structures that can be used to endow an incoming radiation beam with a patterned 
cross-section, corresponding to a pattern that is to be created in a target portion of a 
substrate; the term "light valve" has also been used in this context. Generally, the 
pattern will correspond to a particular functional layer in a device being created in the 
target portion, such as an integrated circuit or other device. 

One example of such a device is a mask, which is generally held by a (movable) 
mask table. The concept of a mask is well known in lithography, and it includes mask 
types such as binary, alternating phase-shifl, and attenuated phase-shift, as well as 
various hybrid mask types. Placement of such a mask in the projection beam causes 
selective transmission (in the case of a transmissive mask) or reflection (in the case of a 
reflective mask) of the radiation impinging on the mask, according to the pattern on the 
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mask. The mask table ensures that the mask can be held at a desired position in the 
incoming projection beam, and that it can be moved relative to the beam if so desired. 

Another example of such a device is a matrix -addressable surface having a 
viscoelastic control layer and a reflective surface. The basic principle behind such an 
apparatus is that (for example) addressed areas of the reflective surface reflect incident 
light as diffracted light, whereas unaddressed areas reflect incident light as undiffracted 

light. Using an appropriate filter, the said undiffracted light can be filtered out of the 
reflected beam, leaving only the diffracted light behind; in this manner, the beam 
becomes patterned according to the addressing pattern of the matrix-addressable surface. 
An alternative embodiment of a programmable mirror array employs a matrix 
arrangement of tiny mirrors, each of which can be individually tilted about an axis by 
applying a suitable localized electric field, or by employing piezoelectric actuation 
means. Once again, the mirrors are matrix-addressable, such that addressed mirrors will 
reflect an incoming radiation beam in a different direction to unaddressed mirrors; in 
this manner, the reflected beam is patterned according to the addressing pattern of the 
matrix-addressable mirrors. The required matrix addressing can be performed using 
suitable electronic means. In both of the situations described hereabove, the patterning 
means can comprise one or more programmable mirror arrays. More information on 
mirror arrays as here referred, to can be gleaned, for example, from United States Patents 
US 5,296,891 and US 5,523,193, and PCT patent applications WO 98/38597 and WO 
98/33096, which are incorporated herein by reference. In the case of a programmable 
mirror array, the said support structure may be embodied as a frame or table, for 
example, which may be Fixed or movable as required. 

Another example is a programmable LCD array, in which case the support 
structure can again be a frame or table, for example. An example of such a construction 
is given in United States Patent US 5,229,872, which is incorporated herein by 
reference. 

For purposes of simplicity, the rest of this text may, at certain locations, 
specifically direct itself to examples involving a mask; however, the general principles 
discussed in such instances should be seen in the broader context of the patterning 
means as hereabove set forth. 
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The term projection system encompasses various types of projection systems. 
Though "lens" in a layperson's understanding usually connotes only refractive optics, 
herein this term is used broadly to include catoptric and catadioptric systems, for 
example. The illumination system may also include elements operating according to any 
of these principles for directing, shaping or controlling the projection beam, and such 
elements may also be referred to below, collectively or singularly, as a "lens". 

Additionally, the term "wafer table" may be used without implying that the 
substrate receiving the image is a silicon wafer, but may rather indicate a stage suited 
for support of any substrate to be processed by the lithography apparatus. 

Lithographic projection apparatus can be used, for example, in the manufacture 
of integrated circuits (ICs). In such a case, the patterning means may generate a circuit 
pattern corresponding to an individual layer of the IC, and this partem can be imaged 
onto a target portion (comprising one or more dies) on a substrate (silicon wafer) that 
has been coated with a layer of radiation-sensitive material (resist). In genera], a single 
wafer will contain a network of adjacent target portions that are successively irradiated 
via the projection system, one at a time. In current apparatus, employing patterning by a 
mask on a mask table, a distinction can be made between two different types of 
machine. In one type of lithographic projection apparatus, each target portion is 
irradiated by exposing the entire mask pattern onto the target portion at once; such an 
apparatus is commonly referred to as a wafer stepper. In an alternative apparatus — 
commonly referred to as a step-and-scan apparatus — each target portion is irradiated 
by progressively scanning the mask pattern under the projection beam in a given 
reference direction (the "scanning" direction) while synchronously scanning the 
substrate table parallel or anti-parallel to this direction. Since, in general, the projection 
system will have a magnification factor M (generally < 1). the speed V at which, the 
substrate table is scanned will be a factor M times that at which the mask table is 
scanned. More information with regard to lithographic devices as here described can be 
gleaned, for example, from US 6,046,792, incorporated herein by reference. 

In a manufacturing process using a lithographic projection apparatus, a pattern 
(e.g. in a mask) is imaged onto a substrate that is at least partially covered by a layer of 
radiation-sensitive material (resist). Prior to this imaging step, the substrate may 
undergo various procedures, such as priming, resist coating and a soft bake. Alter 
exposure, the substrate may be subjected to other procedures, such as a post-exposure 
bake (PEB), development, a hard bake and measurement'inspection of the imaged 
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features. This array of procedures is used as a basis to pattern an individual layer of a 
device, e.g. an IC. Such a patterned layer may then undergo various processes such as 
etching, ion-implantation (doping), metallization, oxidation, chemo-mechanical 
polishing, etc., all intended to finish off an individual layer. If several layers are 
required, then the whole procedure, or a variant thereof, will have to be repeated for 
each new layer. Eventually, an array of devices will be present on the substrate (wafer). 
These devices are then separated from one another by a technique such as dicing or 
sawing, whence the individual devices can be mounted on a carrier, connected to pins, 
etc. Further information regarding such processes can be obtained, for example, from 
the book "Microchip Fabrication: A Practical Guide to Semiconductor Processing", 
Third Edition, by Peter van Zant, McGraw Hill Publishing Co., 1997, ISBN 0-07- 
067250-4, incorporated herein by reference. 

For the sake of simplicity, the projection system may hereinafter be referred to 
as the "lens"; however, this term should be broadly interpreted as encompassing various 
types of projection system, including refractive optics, reflective optics, and catadioptric 
systems, for example. The radiation system may also include components operating 
according to any of these design types for directing, shaping or controlling the 
projection beam of radiation, and such components may also be referred to below, 
collectively or singularly, as a "lens". Further, the lithographic apparatus may be of a 
type having two or more substrate tables (and/or two or more mask tables). In such 
"multiple stage" devices the additional tables may be used in parallel, or preparatory 
steps may be carried out on one or more tables while one or more other tables are being 
used for exposures. Dual stage lithographic apparatus are described, for example, in US 
5,969,441 and WO 9S/40791, incorporated herein by reference. 

As illumination systems have evolved from producing conventional to annular, 
and on to quadrupole and more complicated illumination configurations, the control 
parameters have concurrently become more numerous. In a conventional illumination 
pattern, a circular area including the optical axis is illuminated, the only adjustment to 
the pattern being to alter the outer radius (o r ). Annular illumination requires the 
definition of an inner radius (aj in order to define the illuminated ring. For niuJtipole 
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patterns, the number of parameters which can be controlled continues to increase. For 
example in a quadnipolc illumination configuration, in addition to the two radii, a pole 
angle a. defines the angle subtended by each pole between the selected inner and outer 
radii. 

Concurrently, mask technology has been evolving as well. Binary intensity 
masks have given way to phase-shifted masks and other advanced designs. While a 
binary mask simply transmits, reflects or blocks imaging radiation at a given point, a 
phase-shifted mask may attenuate some radiation or it may transmit or reflect the light 
after imparting a phase shift, or both. Phase-shifted masks have been used in order to 
image features which are on the order of the imaging radiation's wavelength or smaller, 
since diffraction effects at these resolutions can cause poor contrast and end-of-line 
errors, among other problems. 

The various types of illumination configurations can be used to provide 
improvements in resolution, depth of focus, contrast and other characteristics of the 
printed image. However, each illumination type has certain tradeoffs. For example, 
improved contrast may come at the expense of depth of focus; each type of mask has a 
performance which is dependent on the partem to be imaged as well. 

Conventionally, in order to select the optimum illumination mode for a given 
pattern to be imaged onto a wafer, a series of test wafers has been exposed and 
compared on a hit-or-miss basis. As noted above/modern illumination systems have 
ever increasing numbers of variables which can be manipulated. As the various 
permutations of variable settings increase, the cost of trial and error optimization of 
illumination configurations becomes very large and quantitative methods of selecting 
illumination configurations are needed. 

To provide for the above identified needs and others, the present invention 
provides a method of optimizing an illumination profile for a selected patterning means 
pattern, comprising the steps of: 

defining a transmission cross coefficient function for an optical system including 
an illuminator and the selected patterning means pattern; 
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determining relative relevance to imaging of diffraction orders based on the 
selected pattern; and 

calculating an optimized illumination configuration from the transmission cross 
coefficient function, weighting regions of the illumination configuration based on the 
relative relevance to imaging of the diffraction orders. 

According to another aspect of the present invention there is provided a device 
manufacturing method comprising the steps of: 

(a) providing a substrate that is at least partially covered by a layer of radiation- 
sensitive material; 

(b) providing a projection beam of radiation using an illumination system; 

(c) using patterning means to endow the projection beam with a pattern in its 
cross-section; 

(d) projecting the patterned beam of radiation onto a target portion of the layer of 
radiation-sensitive material, wherein, prior to step (d), the cross-sectional intensity 
distribution in the projection beam produced in step (b) is tailored to the pattern 
employed in step (c) using a method as described above. 

According to another aspect of the present invention there is provided a 
lithographic projection apparatus comprising: 

- an illumination system for providing a projection beam of radiation; 

- a support structure for supporting- patterning means, the patterning means 
serving to pattern the projection beam according to a desired pattern; 

- a substrate table for holding a substrate; 

- a projection system for projecting the patterned beam onto a target portion of 
the substrate, 

wherein the apparatus additionally comprises: 

- calculating means, for defining a transmission cross coefficient function for the 
illuminator and the desired pattern, determining relative relevance to imaging of 
diffraction orders based on the pattern produced by die patterning means, and 
calculating an optimized illumination configuration from the transmission cross 
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coefficient function, weighing regions of the illumination configuration based on the 
relative relevance to imaging of the diffraction orders; and 

- selecting means, for selecting the cross-sectional intensity distribution in the 
projection beam exiting the ulumination system in accordance with the illumination 
configuration calculated by the calculating means. 

According to a yet another aspect of the present invention there is provided a 
method of optimizing a selected mask design comprising: 

identifying critical features of the selected mask design; 

detennining an optimized illumination profile based on diffraction orders of the 
critical features; and 

modifying the selected mask design by use of optical proximity correction 
techniques which are selected to reduce a number of pitches present in the selected mask 
design. 

The present invention further provides computer programs for carrying out the 
methods described above. 

Although specific reference may be made in this text to the use of the apparatus 
according to the invention in the manufacture of ICs, it should be explicitly understood 
that such an apparatus has many other possible applications. For example, it may be 
employed in the manufacture of integrated optical systems, guidance and detection 
patterns for magnetic domain memories, liquid-crystal display panels, thin-film 
magnetic heads, etc. The skilled artisan will appreciate that, in the context of such 
alternative applications, any use of the terms "reticle", "wafer" or "die" in this text 
should be considered as being replaced by the more general terms "mask", "substrate" 
and "target portion", respectively. 

The present invention will be described further below with reference to 
exemplary embodiments and the accompanying drawings, in which: 

The present invention involves first mathematically modeling the imaging of the 
pattern onto the substrate (e.g. from a mask), taking into account the illumination source 
and the pattern details. 
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There are two primary methods for calculating the aerial image for a finite 
illumination source. These methods are Abbe's formulation and Hopkins's formulation. 
In Abbe's formulation, each point source in the illumination configuration produces a 
plane wave incident onto the pattern, and each of these point sources is imaged onto the 
wafer. Since the source is spatially incoherent, the total intensity at the wafer is the 
summation of the intensity produced by each of these point sources. Therefore, in 
Abbe's formulation, the integration over the illumination configuration is performed 
after the integration over the pattern. 

In Hopkins's formulation, the order of integration is switched, i.e., the 
integration over the source is performed first. In Hopkins's formulation, a four 
dimensional transmission cross coefficient (TCC) is defined, and the image intensity is 
the inverse Fourier transform of the TCC. A derivation of the TCC is described, for 
example, in Bom and Wolf, Principles of Optics , 6" 1 Ed., pp. 528-532, herein 
incorporated by reference. 

The TCC is the autocorrelation of the projection pupil multiplied by the 
ilhunination pupil. The TCC is shown in Figure 1 as a set of three overlapping circles. 
From left to right, the first circle represents the illumination pupil J,(a,p) where a and p 
are coordinates of the illumination configuration. For the purposes of the following 
calculations, the radius of J, may be, for example, set to the maximum allowable outer 
a, for the lithography apparatus which will be used for imaging. It is also possible to set 
cr r to 1.0 or larger in order to perform feasibility studies and to determine the benefits of 
a larger ov 

The central circle represents the projection pupil, K(a,P) that is centered at 
(-mA/P,NA, -nA/PyNA). The coordinate systems arc normalized by a factor of A/NA so 
that the radius of K is 1.0. The circle on the right likewise represents the projection 
pupil; however, it is centered at (pA/P.NA, qX/PyNA). In these last two expressions, 
m,n,p, and q correspond to discrete diffraction orders and it becomes clear that the TCC 
is a four dimensional (4-D) equation as described above. The diffraction orders in the x- 
direction are represented by m and p and the diffraction orders in the y-direction are 
represented by n and q. Though for puiposes of this description, x and y coordinates are 
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used, one skilled in the art understands that alternate coordinate systems could be used 
with appropriate changes of coordinate systems in the following equations. 

The TCC for a 4-D discrete point (m,n,p,q) is the integral of the shaded area 
where all three circles overlap. Since the structure is assumed to be periodic, the 
Fourier transform of the pattern is discrete and the TCC is discrete. For a continuous 
pattern image, the pitch can be increased until an adjacent feature has no influence on 
the Fourier transform of the partem of interest. The TCC in Figure 1 is described 
mathematically in Equation 1 . 

TCC(m,n,p.q) = 

The TCC may be expanded to include the effects of the pattern by defining 
diffraction order cross coefficients (DOCC). The DOCC are defined in Equation 2 
which is derived from the multiplication of the TCC by the Fourier transform 
coefficients of the partem. 

DOC<Xm,n,p,q) = T(m, n y*{-p,-q)rCC(m,n,p,q) Eqn. 2 

Further, the radiation intensity at the wafer may be calculated by the inverse 
Fourier transform of the DOCC, as shown in Eqn. 3. 

H^.y)-ZZXZ e L " ] * [ * ] DOCC{n,.n.p.q) Eqn. 3 

m * p q 

The projection optical system acts in part as a low pass filter, which reduces the 
diffraction orders so that only a few of the diffraction orders are important to the 
calculated image intensity. As a result, the TCC is a band limited function. The 
maximum necessary x and y orders can be calculated according to Equations 4 and 5 
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respectively. In each case, both the negative and positive orders are necessary, for 
example m extends from negative to positive nu.^, < m < +m mm>x ). Since 
both negative and positive orders are needed, the size of the TCC is 21x^+1 by 2n imx +l 
by ^max" 1 "! by 2^+1. Fortunately, however, because the TCC is band limited, it is not 
necessary to calculate all the pattern diffraction orders. Like in the TCC, only partem 
diffraction orders -m^^ m < +m m „ in the x direction and orders -n^S n < +11^ in the 
y direction are needed. 



Substituting Equations 1 and 2 into Equation 3 gives Equation 6 for the radiation 
intensity at the wafer. By switching the order of integration, as shown in Equation 7, 
that is, by using Abbe's formulation rather than Hopkins's, the portions of the 
illumination pupil which are most influential on imaging may be determined. Note that 
each of Equations 6 and 7 extends across two lines. 




Eqn. 4 




Eqn. 5 




Eqn. 6 
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Since a and JJ represent illumination pupil coordinates, a new function, J,,,,, may 
be defined. The new function J opl indicates which part of the illumination configuration 
(ct,(3) is being used for a given diffraction order (m,n,p,q) and is expressed in Equation 
8. From Equation 8, the image intensity can be calculated by multiplying it by the 
inverse Fourier coefficient (e ta ) and summing over all 6 variables (m,n,p,q 1 a,p) as 
shown in Equation 9. 



Eqn. 8 



I(x,y)= tf^PZIZZ^"^ Eqn. 9 



As will be appreciated, J opl is a six dimensional function and it is therefore 
difficult to apply it to the illumination configuration. In order to best determine which 
portions of the illumination configuration are significant to image formation, it is 
desirable to eliminate a few of the six variables. 

The aerial image intensity, l(x,y), is found by taking an inverse transform over 
m+p and n+q. When m+p=n+q=0, there is no modulation in the aerial image intensity. 
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Since one of the goals of illumination optimization is to eliminate parts of the 
illumination configuration that have little or no influence on modulation, those portions 
of the illumination configuration for which m+p = n+q = 0 may he ehminatcd. In order 
to eliminate these parts and to better visualize the illumination configuration portions 
significant to image formation, a transformation of variables will eliminate two of the 
variables in the six dimensional J opt function (four diffraction orders) and convert it into 
a four dimensional function (two diffraction orders). The four dimensional function is 
called J m<tD . By substituting Equations 10 and 1 1 into Equation 9 for I(x,y), Equation 
1 2 may be derived. 

rj = m+p=> p = r)-m Eqn. 10 

£ = n + <j=>ej = f-n Eq n H 



Eqn. 12 

In Equation 12, J,,^ can be seen to be the summation of J,,,, over m and n after 
the transformation of variables according to Equations 1 0 and 1 1 . By further 
substituting Equation 8 into Equation 12, J,,^ may be expressed as in Equation 13 and 
the intensity, I(x.y), can be written as a function of J op ,. 2D as in Equation 14. 



J a p,-lD[a,P,ri,S) = J,{a,fi) ^g™ 'jT T{r»,n)r'[-{r] - m\-{£ - „)] Eqn. 13 

'•-/,».»- /.no. 
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opt-tD 



Eqn. 14 



The function, J opl . 2D , when evaluated, shows those portions of the illumination 
configuration that are important for each diffraction order. Since J^,.^ is weighted by 
each diffraction order, T(m^i), large diffraction orders will have a greater influence on 
the aerial image. 

A starting point for the best illumination configuration for a particular pattern 
may be denoted J tot and is found by summing J,,,,,.^ over t\ and ^ and subtracting 
Jop«d(<*»P >Tl=0,c>=0) as shown in Equation 15. In Equation 15, when n=0 and c=0 there 
is no modulation in the aerial image and the J opt . J1) (a.,p\T 1 =0 > c;=0) component represents 
zero order or DC light. Points in the illumination that do not contribute to imaging 
increase the overall amount of DC light. Since the increased DC light causes no 
modulation it is not of great benefit, and moreover it can result in a reduction in depth of 
focus. 

Thus, an illumination configuration in accordance with J lol minimizes the amount 
of DC light and results in an improved process window. The equation J m can be used to 
show which parts of the illuminator are more significant (or less significant) to image 
formation. 



Eqn. 15 



Since the illumination configuration and pattern are coupled, optical proximity 
correction (OPC) changes influence the diffraction orders, which therefore influences 
J lol . Consequently, modifications to the starting illumination configuration, J to „ and the 
pattern should be performed a few times using iterations of processing with an OPC 
engine and an illumination engine, as will be understood by one skilled in the art. 
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Furthermore, the pattern and illumination configuration also need to be adjusted to 
optimize a particular imaging criteria (depth of focus (DOF), end of line (EOL), 
aberration sensitivity, etc.) which may be performed with optimization software. 
However, since the pattern as a whole, rather than the OPC features, has the largest 
impact on the optimal ulumination configuration, J** is the best initial illumination 
configuration to lead to the fastest convergence for optimizing iterations over the 
illumination configuration and pattern. 

The starting illumination configuration, J B1 , may be represented by a gray scale 
illumination configuration having continuous values of intensity over a 0 to 1 range. It 
is possible to create such a gray scale illumination configuration with a diffiactive 
optical element (DOE) or by using a quartz plate with dithered chromium plating. If a 
gray scale illumination configuration is not possible or preferred, the illuminator profile 
can be forced to be only 0 and 1 by applying a threshold to the gray scale, in which 
values above the threshold are rounded up to 1 and values below the threshold are 
rounded down to 0. An arbitrary threshold can be applied, or an optima] threshold may 
be found through simulating the process window, or by repeated test runs. 

Example 1 : The technique for calculating J^, outlined above was applied to a 
brick wall isolation pattern. A ISOnm pattern was shrunk to 130nm and 1 lOnm design 
rules and imaged with a step and scan lithography system having a numerical aperture 
(NA) of 0.8. The isolation pattern for the 130nm design rule is shown in Figure 2. 

The magnitudes of the diffraction orders of this mask feature are plotted in 
Figure 3. In Figure 3, the largest order is the (0,0) order or the DC background light. 
The orders that contribute the most to imaging are the (±2,0) orders and represent the 
vertical bricks in the brick wall pattern. The other significant order is the (±1,±1) which 
represents die clear areas and defines the end of the isolation pattern. The higher orders 
also help to define two dimensional structures such as the end of each line. Since the 
diffraction orders are not constant, the orders change the weighting coefficients in the 
DOCC, which implies that the mask pattern influences the illumination strategy. 

The diffraction order coefficients T(m,n) in Figure 3 can be substituted into 
Equation 13 to calculate and are plotted in Figure 4. As can be seen from Figure 
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4, the largest contribution to J,^ is the (n^c^O) order. The (0,0) order does not 
contribute to imaging and decreases the DOF. As Equation 15 shows, this (0,0) order 
can be subtracted from the total illumination, J^,. Not considering the (0,0) order, the 
largest contribution is the (r|=±2,£=0) diffraction order, which represents the formation 
of the isolation lines along the x-direction. Another component that is large and defines 
the end of the isolation lines is the (n=± 1,£=^±1) diffraction order. Although the (0,±2) 
diffraction order is rather small, higher orders combine in the t|=0 and £=±2 region of 
the lens. These regions also help to define the end of line. The DOCC approach shows 
how the illumination pupil is sampled for improved image formation and is an effective 
method for understanding imaging of the brick wall isolation pattern. 

Using Equation 15, the iUumination pupil of the brick wall pattern for the 130nm 
design rule can be calculated and is shown in Figure 5. Figure 5 shows that the most 
significant areas for image formation are the outer portions of the illumination 
configuration along the x-axis. These outer portions form an elliptical dipole. In 
addition to these elliptical dipole elements, the center of the illumination pupil has a 
large contribution to image formulation. As noted above, the illumination pupil can be 
implemented in gray scale or binary illumination profiles. 

Depending on the apparatus being used, gray scale illumination may be possible. 
By gray scale illumination is meant controllable illumination intensity for which a 
normalized level ranging from 0 to 1 may be selected for at least given portions of the 
illumination configuration. For example, such control over illumination intensity may 
be produced by use of a diffractive optical element (DOE) in the illumination system. 
In this case, for example, the illumination configuration can be implemented as shown 
in Figure 5. However, some of the local spikes which are calculated in theory and seen 
in Figure 5 will be removed after low pass filtering of the illumination information as a 
result of the projection optics, as discussed above. Therefore, when designing an 
illumination configuration, spikes which will be filtered should be ignored. 

If a binary illumination configuration is used, i.e., only binary values for 
intensity of the illuminator are allowed (0 or 1), a threshold value should be chosen as a 
basis for assigning values of 0 or 1 to each point on the illumination configuration. For 
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example, if a threshold of 0.8 is chosen, illuminator intensity values above 0.8 are 
rounded up to 1 and values below 0.8 are rounded down to 0. Other threshold values 
may be applied as desired. 

Example 2: Using the gray scale to binary approach, a binary illumination 
configuration for the same brick wall isolation partem was designed assuming a 
maximum outer c of 0.88 and is shown in Figure 6. 

The performance of the optimized illumination configuration in Figure 6 was 
then simulated for binary mask on a step and scan photolithography apparatus having 
NA=0.8 and k=248nm and compared to the simulated performance of annular 
ulurrtination. In the simulation, the vector (thin-film) imaging resist model was used 
since the numerical aperture was above 0.7. In this model, the resist is 400nm thickness 
of a type having a refractive index n=I.76-j0.01 16), over 66nm of another type having 
n=1.45-j0.3 on top of a polysilicon material having n=l .577-j3.588. The results with 
the annular illumination (^=0.58 and of a ou( =0.88) and with the optimized iUuminator 
(a^ =0.88) are shown in Figures 7 and 8, respectively. In both Figures 7 and 8, cross 
section results in the middle of the isolation region and top down simulation results are 
shown. In the Figures, the Bossung plot B from aerial image threshold is calculated by 
averaging the intensity through the resist, and the resultant linewidth, lw, is plotted 
versus focus, f, for a threshold intensity. This technique tends to over predict the DOF 
as thickness loss and resist profile slope is not taken into consideration. A resist model 
that at least calculates thickness loss is probably necessary. In each of the figures, the 
top down results are plotted as solid, curved lines at the best threshold (best dose) as 
calculated by the Bossung plot. These simulated threshold images are compared to 
actual mask data shown in dashed, straight lines. 

Simulation results for the 130nm design rule brick wall isolation pattern are 
plotted in Figure 7 for a binary mask feature with NA of 0.8 using annular illumination 
(o iu =0.58 and of o o „,=0.88). This annular setting has approximately 0.4um of DOF from 
-0.4j.im to O.Oum focus. The contrast of the resist is low through focus, and can be 
imaged with a low contrast resist. However, at this low intensity conn-ast, the mask 
enor enhancement factor (MEEF) is large and the exposure latitude (EL) is small. The 



(43) 



^12002-3 3483 6 



top down images in Figure 7 also show that there is approximately 20nm of end of line 
(EOL) shortening, which can be fixed by extending the line slightly for the 130nm 
design rule. However, as the design rule continues to shrink, extending the line is no 
longer feasible as the extended line may conflict with other features. Therefore, it is 
desirable to fix the EOL with the illumination. 

In Figure 8, simulation results for the 130nm design rule brick wall isolation 
pattern are plotted for a binary mask feature with NA of 0.8 and using the optimized 
binary illumination configuration in FIG. 6. The optimal fllumination configuration has 
approximately 0.6um DOF from -0.45 pm to +0.15|im focus. In comparing the cross 
section images in Figure 8 to those in Figure 7, the optimized illumination configuration 
has a larger contrast through focus as compared to annular illumination. This larger 
contrast implies that the MEEF for the optimized illumination configuration is lower 
compared to annular illumination and that the exposure latitude for the optimized 
illumination configuration is higher. Another benefit of the optimized illumination 
configuration is the improved line end performance as compared to annular 
illumination. The top down images in Figure 8 show that the optimized illumination 
configuration is capable of maintaining the EOL without extending the line on the 
pattern, which is advantageous for more aggressive design rule shrinks. 

Example 3: The results in Figures 7 and 8 for binary mask (BIM) were 
compared to simulation results for chromeless mask (CLM). A chromeless brick wall 
isolation pattern was designed from experimental results of software simulation in a 
manner known to those skilled in the art. The chromeless technology requires (0,0) 
order light so as to fully benefit from the DOF improvement produced by off axis 
illumination. Experimental results from the simulation confirm the need for (0,0) order 
light for which purpose die isolation layer should be dithered or half toned. The half 
tone pitch may be chosen such that the first order in the dithered direction does not fall 
into the projection pupil. In the example, the lines were dithered in the vertical direction 
with pitch less than Xy[NA(l+a otI1 )]. The dithering duty cycle however should be tuned 
to optimize the amount of (0,0) order light for best DOF and pattern fidelity. In the 
simulation results for CLM the half tone pitch was 155nm with 50% duty cycle (77.5 run 
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chrome islands). This pitch substantially prevented the (0+1) orders from entering the 
projection pupil; however, this duty cycle should be tuned to maximum DOF with 
computer aided design tools. 

Example 4: Simulation results for the 130nm design rule layer were plotted for 
a CLM with 155nm halftone pitch and 50% duty cycle. The CLM was exposed on a 
X=248nm apparatus with NA of 0.8 and annular illuminauon (c in =0.58 and o out =0.88). 
The CLM with this annular setting had 0.5pm DOF (-0.4pm focus to +0.1pm focus). 
The CLM with annular illumination had greater DOF and better contrast through focus 
as compared to BIM with' annular illumination. This indicates that the CLM performed 
better than the BIM mask. The top down simulation results indicated that the EOL 
performance with CLM is theoretically better than the EOL performance with BIM and 
that the CLM was also able to better define the contact hole landing area as compared to 
BIM. 

Example 5: Simulation results for 130nm brick wall isolation pattern isolation 
layer were plotted for X-248nm apparatus with NA of 0.8 and the optimized elliptical 
dipole shown in Figure 6. These results were simulated with a reticle identical to the 
CLM reticle used in the preceding example, which has a 155nm halftone pitch and 50% 
duty cycle. The CLM exposed with this optimized illumination configuration had 
0.7um DOF (-0.5pm to +0.2pm), an improvement of 40%. The Bossung plots indicated 
that the isofocal intensity was around 0.21. A model based OPC approach could 
additionally be applied in order to tune the reticle to size at the correct linewidth, 
providing further improvements in performance. The linewidth may be corrected, for 
example, by biasing and by modifying the half tone duty cycle. The top down 
simulation results indicated that CLM was able to define the contact landing region and 
was able to maintain CD uniformity. Necking and other linewidth inconsistencies were 
reduced with this elliptical illumination configuration. Furthermore, the CLM reticle 
could be biased to improve the DOF, and consequently the EOL perfonnance should 
improve. Furthermore, model based OPC should be able to correct the EOL further. 

Example 6: Using the mask pattern of Figure 2 for a 1 lOnm design rule 
isolation layer, an optimized illumination configuration was generated with Equations 
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13 and 15. In order to visualize sampling of the illumination pupil, J^,.^ is plotted in 
Figure 9, showing x orders (r)=m+p) horizontally and y orders (£=n+q) vertically. As in 
Figure 4 for the I30nm design rule, the largest contribution to 1 lOnm design rule in 
Figure 1 1 is the (rj=0, £,=0) order. This (0,0) order light is detrimental to the DOF and is 
eliminated in J ttt as indicated in Equation 15. Figure 9 also shows that the (±1 ,±1) 
orders are the largest contributors to the illumination configuration optimization rather 
than the (±2,0) order. This is due to the fact that 1 1 Onm design rule is too aggressive for 
the 24Snm apparatus with NA=r=0.8, as slightly higher NA is preferred to achieve this 
resolution. The orders that contribute most to defining the isolation linewidth are the 
(±2,0) orders. The (±2,0) orders, however, are at the far edge of the illumination 
configuration (0.8<ct<1.0), which indicates that a of 1 may provide improvement to 
implement the 1 lOnm design rule at this wavelength. 

Using Equation 15 and the results in Figure 9, the optimized illurriination 
configuration for the 1 1 Onm brick wall isolation layer is shown in Figure 10. Figure 10 
shows that the illumination configuration areas that contribute most to image formation 
are a small portion in the center and far edges of the illumination configuration. One 
possible implementation of this illumination configuration is plotted in Figure 11a. In 
order to use 248nm apparatus to print more aggressive design rules and push the limit of 
the projection numerical aperture, an illumination configuration with a of 1 .0, as shown 
in Figure 1 lb, and with small sectors (ct ring width of 0.2) may be used. 

An implementation of the present invention includes selection of cells or 
particular gates that are critical. These critical features are then processed to determine 
J 10t as described above. In Section 1 , it was shown that the illumination configuration is 
pattern dependent. Therefore, if the pitch does not differ significantly for the critical 
features, it is possible to create a single illumination configuration which optimizes the 
process window for all the critical features. In Figure 12, an example of a circuit with 
critical gates g„ g 3 , g 3 and a critical cell cc is shown. The diffraction orders of these 
tagged critical features can be calculated, and by using the theory already described, the 
optimized illumination configuration can be calculated. After calculating the optimized 
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illumination configuration the process window can be calculated and compared to the 
process window with other iUumination configurations. 

Another method of optimizing iUumination/pattern interaction is to modify the 
pattern design with scattering bars. Scattering bars discretize the pitch from a semi- 
continuous function for an ASIC or logic design. After placing scattering bars, there are 
fewer pitches. This can be demonstrated in simulation software by placing the 
scattering bars with an edge-to-edgc separation of 0.61X/NA. In Figure 13, the design 
in Figure 1 2 has been modified by adding a plurality of scattering bars. The 
illumination configuration can then be optimized for the modified design. The process 
window performance of an illumination configuration optimized for a design with 
scattering bars can then be compared to the process window of an illumination 
configuration optimized without scattering bars. Since a design with scattering bars 
discretizes the pitch, the combination of scattering bars with optimized off axis 
illumination (OAI) will have the largest possible DOF process window. 

Another concept for implementing illumination configuration optimization is 
through the placing of scattering bars based on space width (SW) considerations. A 
scattering bar is placed through rule based OPC which rules may be denned by the 
space width. With simulation software, it should be possible to calculate the probability 
density function (pdf) of the space width without scattering bars and with scattering 
bars. The illumination can than be optimized by considering the pdf by modifying 
as shown in Equation 16. Assuming that the vertical lines and horizontal lines are 
infinite, it is also possible to calculate the diffraction orders T(m,n). In Equation 17, the 
diffraction orders are calculated as a function of m and n where w is the line width, x is 
the intensity transmittance of the reticle, and P,=S W,+w and P y =SW y +w are the pitches 
in the x and y directions, respectively. 
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Equation 17 is a matrix of four equations for which, in order of presentation, 
m=n=0; m=0, a£0; m?£0, n=0; and m^Q, n^O. 



Calculating the optimal illumination configuration with the pdf can present some 
problems as it implies that some pitches are not as important as others. If all the gates in 
the pdf are considered to be critical, the pdf should be modified by a weighting factor. 
This weighting factor is a function of pitch called wf(Px). With this weighting factor, 
all the critical pitches should be treated the same such that wf(Px>pdf(Px)=l . This 
weighting factor should be added to Equation 16 by replacing pdfTPx) in Equation 16 
with wf(Px)»pdf(Px). In the case that all of the pitches are critical, the weighting factors 
will not help to resolve the optimization and it is difficult to create an optimized 
illumination configuration without modifying the (pattern) design. 

One solution to this problem is to modify the design by adding scattering bars as 
discussed above. Scattering bars help to reduce pitches for isolated features. Once 
scattering bars are added to die design, the previously isolated features tend to act like 
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dense features. Thus, scattering bars discretize the pitch from a continuous pdf to a 
more discrete pdf. Figure 14 is an example pdf for a logic pattern with features oriented 
in the y-direction (i.e. "vertical" direction) in which scattering bars have and have not 
been applied. Figure 14 shows vertical gate space widths (um) on the x (horizontal) 
axis. For the unmodified design, D, without scattering bars, there are three discrete 
humps in the pdf at space widths of 0.2, 0.6 and 1.5nm. After the placement of 
scattering bars, D+SB, the number of pitches has been reduced such that most of the 
space widths are at dense pitches of 0.2um. With this change to the pdf, it is more 
likely that an illumination configuration can be optimized. 

The total illumination configuration for a design with both horizontal (x- 
direction) and vertical features is the sum of horizontal and vertical illumination 
configurations. If the illumination configuration is concentrated at ct„ for the vertical 
features and is concentrated at <j cy for the horizontal features, the optimal illurnination 
configuration will be a "conventional" quadrupole illumination configuration provided 
that v5o-„ £ l and that VSo-g, <: l . Otherwise, this type of analysis will result in a four 
pole illurnination configuration that has been rotated 45° 

The illumination technique presented herein maybe extended to account for 
aberrations. Including aberrations allows an operator to determine which part of the 
illumination configuration couples to the aberration. The amount of coupling is directly 
related to sensitivity of the image intensity to the aberration. By understanding this 
coupling, it may be possible to modify the illumination configuration to minimize the 
aberration sensitivity of a design. 

The projection pupil, K(a,3), for scalar imaging contains the obliquity factor, 
defocus, and the exponential of the wavefront represented by the Zemike polynomials. 
This scalar imaging pupil is shown in Equation 18. This pupil can further be divided 
into two parts, the imaberrated pupil K^a.P) and the aberrated pupil (the exponential of 
the wavefront); these two parts are multiplied together as shown in Equation 19. 
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Eqn. 18 



K{a, fl) = K 0 (a,fl)eJ-i2Z w{a,fi)\ 

unaberrateei y — 3 



aberrations 



Eqn. 19 



where 



obiltqutry- factor 



defocus 



Eqn. 20 
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Eqn. 21 



* 2 .t 3 



e =I+x+ — + — + ... = ] +. .x 
2! 3! 



Eqn. 22 



From Equation 22, the wavefront can be written as a linear approximation, 
which is shown in Equation 23. By substituting Equation 23 into Equation 22, the 
linear approximation for the projection pupil, K(a,(J), can be calculated with Equation 
24. 



e«p[-/^»K(«.^] a i-<^lK«./?)=l-/^g^jr K (« 1 >») Eqn. 23 
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K{a,fi) = K 0 { a ,j3^-i~J:z v R v (a,fi)^ Eqn. 24 

Since the TCC is a function of the projection pupil, K(a,P), the linear 
approximation to the pupil in Equation 24 implies that the TCC can be represented by a 
linear approximation. This is accomplished by substituting Equation 24 into Equation 
1 , which results in Equation 25 . Once again by neglecting the terms of power 2 or 
greater, the TCC in Equation 25 can be simplified as shown in Equation 26. 

The wavefront, W(a,(3), is most often denoted by the summation of Zernike 
fringe polynomials as shown in Equation 2 1 . Using the linear theory of aberrations, the 
exponential, e x , can be represented by a Taylor series expansion. The Taylor series 
expansion is valid for small x, and previous work has shown good agreement for aerial 
images when Z v is less than 0.04k. The Taylor series expansion for e" is shown in 
Equation 22. In Equation 22, terms of power 2 or greater have been dropped, which is 
valid provided that Z v is less than 0.04 (0.04 2 =0.0016 and is negligible). 




Eqn. 25 
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Eqn. 26 
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By defining the unaberrated TCC, TCCoOn.n.p.q), and the aberrated TCC, 
TCC v (m,n,p,q), in Equations 27 and 28, respectively, the TCC can be represented by a 
linear function of TCC 0 and TCC V as shown in Equation 29. 

Eqn. 27 



TCC v ( m ,n, P ,g) = -i^- JJ a + j£L,fi+J!±-] . Eqn. 28 



TCC(m, n,p.q) = TCC 0 {m,n, p, g)+ £ Z„ ^CC„ (m, »,/>,?)+ TCC^- p ,-g ,-m,-/i)] 

v-5 

Eqn. 29 

Since the TCC can be constructed as a linear approximation as shown in 
Equation 29, J opt can also be written as a linear approximation. The linear 
approximation to is derived in Equation 30 by using Equation 8 for J,^ and by 
following the methodology for the linear approximation of the TCC as outlined in 
Equations 18 through 29. 
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Eqn. 30 



Equation 30 for J opl can then be divided into a summation of the unaberrated 
with the aberrated J^,* as shown in Equation 33. The definition of J^^ and J oplV are 
shown in Equations 31 and 32, respectively. 
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Eqn. 31 
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Eqn. 32 
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Equation 32 describes the portion of the illumination configuration which 
couples to a particular aberration. The amount of coupling affects the image intensity 
and helps to provide an understanding of aberration sensitivity to illumination. By 
combining Equations 31 and 32, J opt can be written as a linear approximation. 

J opl (a, 0, m, n, p. q)s J opl0 (a, 0. m,n, p, q)+ £ 2 V [/^(cr, 0,m, n, p,q)+ j' opni {a, 0,-p ,-q,-m,-«)] 



v=5 



Eqn. 33 



In another aspect of the present invention, weighting factors can be introduced to 
maximize or to minimize a response to a particular metric including, for example, depth 
of focus (DOF), image log slope (ELS), image slope (IS), or aberration sensitivity. The 
optimal J*,, of Equation 15 can be modified to include these weighting factors as shown 
in Equation 34. 

a 11 p 7 



In general, photoresists react in proportion to the logarithm of the intensity of the 
light impinging thereon. As the intensity, and therefore the logarithm of the intensity, 
increases, the feature will be printed into the resist with better fidelity (i.e. improved 
resist profile and improved process window). Therefore, it is desirable to maximize the 
logarithmic change in intensity (ILS). The 1LS is defined in Equation 35. 
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„ „ dlnl 1 3/ . 

ILScc—— = -— Eqn. 35 

or i etc 



Since the derivative of the intensity changes faster than the inverse of intensity, 
Equation 35 will increase more by increasing the derivative of the intensity. The 
intensity can be calculated from Equation 3 and the derivative of the intensity with 
respect to x is defined in Equation 36. The derivative with respect to x results in the 
weighting function, w x , as shown in Equation 37. Likewise, a weighting function w y 
can be defined with respect to y as shown in Equation 38. 

- ZZZP^O" 'beet*-.**) 



w. 

Eqn. 36 

vv r = i — (m 4- /?) Eqn. 37 

w v =i'^(« + g) Eqn. 38 



Since pattern features and intensity features are two dimensional, the norm of the 
gradient may be used to indicate the change in intensity with respect to position. The 
norm of the intensity gradient is denned in Equation 39. This allows us to define a 
weighting function to calculate J m in Equation 34. The weighting function to maximize 
the image log slope is defined by Equation 40. 



w IUi {m,n,p,q) = 2x. 
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Eqn. 39 
Eqn. 40 
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Equation 40 shows that when m+p=0 and n+q=0, the weighting function 
becomes 0. When m+p=0 and n+q=0, these orders contribute nothing to image 
modulation and reflect DC contributions to the image. Furthermore, w^ increases as 
m+p and n+q increase. This implies that higher order diffraction order terms are more 
highly weighted and contribute more to ILS. 

In addition to maximizing ILS, the depth of focus of the process will increase if 
the ILS is improved such that the intensity response to focus is minimized. The focus is 
accounted for by the pupil K(a,p). The pupil, K(a,|5), is shown in Equation 41, where 
focus is denoted z. Equation 41 can be divided into two terms, terms dependant on z 
(the defocus term) and terms independent of z (the non-defocus term), as shown in 
Equation 42. 



i-(a 2 +P 2 ) J ^[" T^ " p J«4- t (a ' W J qn ' 41 
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~ defocus ed 



~K m/ (a,/!)KAa,fi) Eqn. 42 



The variation in intensity due to focus, z, can be minimized by setting the 
derivative of intensity with respect to z to zero. By substituting Equation 42 into 
Equations 1 through 3, a cost function f(cc,|3,z) can be defined as shown in Equation 43 
which is the cost function of the intensity imaging terms which depend on z. 



f{a,P,z)=K a 
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Eqn. 43 



The cost function, f(a,P,z), is in turn minimized when g(a,p\ > m,n,p,q) equals 
zero (see Equation 44, below). In Equation 44 the phase terms have been removed as 
the derivative with respect to z equals zero only when the magnitude terms equal zero. 
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When g(a,p , ,m,n,p,q) is zero, the areas (a,f}) of the pupil for a given order (m,n,p,q) are 
minimally sensitive to focus. These are the most desired areas of the pupil for 
constructing the illumination configuration. A weighting function, w foiM (a,(},m,n,p,q), 
is defined in Equation 45. This weighting function is equal to 1 for areas least sensitive 
to focus and is equal 0 for areas most sensitive to focus. A new weighting function 
which maximizes ILS through focus can then be defined by Equation 46 and used to 
modify the illumination configuration. 



g{a, J9, n, n, p, q) = — f{a, fi, z) 

02 



Eqn. 44 



w /«,« ( a > P> w < "» P > *) = 1 ~ k( a ' A w» »• P> l\ 



Eqn. 45 



w(a, fi,m,n,p,q)'=w lls (m,n, p, q)w /oatt (a, fi, m, n,p,q) 



Eqn. 46 



The above methodology allows the sensitivity of the intensity to be mhiimized to 
the influence of focus, an aberration. Since the impact of focus on intensity can be 
minimized, the impact of intensity can be minimized to a specific aberration. This is 
desirable for certain patterns which demonstrate a high sensitivity to a particular 
aberration. The projection pupil in Equation 19 can be written as an unaberrated term, 
Ko(a,(3), multiplied by an aberrated term, K,(a.3) as shown in Equation 47. 



Eqn. 47 



The sensitivity of the intensity to a particular aberration Z, can be minimized by 
setting the derivative of the intensity widi respect to Z ; to zero. By substituting 
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Equation 47 into Equations 1 through 3 and taking the derivative of the intensity, the 
aberration sensitivity is minimized when h(a,p,m,n,p,q) in Equation 48 is equal to zero. 
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Eqn. 48 



h{a,fi, m,n,p,q) = 
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= 0 
Eqn. 49 



Equation 48 may be simplified and written as Equation 49. A weighting 
function w a6 (a,p^n,n,p,q) is defined in Equation 50, which equates to 1 for areas (ct,p) 
of the pupil that are least sensitive to Zj and to 0 for areas most sensitive to Z;. 



Eqn. 50 



The weighting function to minimize US sensitivity to a particular aberration, Zj, 
can then be defined in Equation 51 . Furthermore, a weighting function to minimize ILS 
sensitivity to a particular aberration, Z,, and to maximize ILS through focus can also be 
defined in Equation 52. Either of these equations can be substituted into Equation 34 to 
calculate the illuminator with the optimal response to a given metric. 



w(a,fi,m,n,p,q) = w NILS [m,n,p,q)w all (a, fi,m, n, p, q) 



Eqn. 51 
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w(a, p, m, n, p, q) = WmL5 (m, n, p.q)*^ {a, P, m, n, p, q )w atl (a, p, m, n, p,q) 

Eqn. 52 

Figure 15 is a schematic representation of an example of a lithography apparatus 
for use according to the present invention. The apparatus includes a radiation system. 
The radiation system is made up of a lamp LA (which may be an excimer laser, for 
example) and an illumination system which may comprise beam shaping optics EX, an 
integrator IN, and a condenser lens CO, for example. The radiation system supplies a 
projection beam PB of radiation. For example, the radiation system may provide 
ultraviolet, deep ultraviolet or extreme ultraviolet radiation. In general, the radiation 
system may also provide soft x-ray or other forms of radiation. 

A first object table, or mask table MT holds a mask MA. The mask MA 
includes a pattern area C which contains the mask pattern to be imaged. The mask table 
MT is movable relative to the projection beam PB so that different portions of the mask 
may be irradiated. Alignment marks M,and M 2 are used for detennining whether the 
mask is properly aligned with the substrate, or wafer, W. 

A projection system PL projects the projection beam PB onto the wafer W. The 
wafer W includes two alignment marks P, and T 2 which are aligned with the marks M, 
and Mj prior to beginning imaging. The wafer W is supported by a substrate table WT 
which is moveable relative to the projection beam for exposing different parts of the 
wafer W; in this way, the mask pattern C may be imaged onto different target portions c 
of the wafer W. An interferometric position monitor IF is used to insure that the wafer 
table WT is in the correct position relative to the position of the mask table MT. 

While the invention has been described in connection with particular 
embodiments, it is to be understood that the invention is not limited to the disclosed 
embodiments, but on the contrary it is intended to cover various modifications and 
equivalent arrangement included within the scope of the claims which follow. 
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4 Brief Description of Bracings 

Figure 1 is a diagram of the transmission cross coefficient function for a 
generalized image forming system; 

Figure 2 is an example of a brick wall isolation pattern microlithographic mask 

feature; 

Figure 3 is a representation of the diffraction orders of the mask feature of 
Figure 2; 

Figure 4 is a map of the calculated optimized four dimensional illumination 
configuration for the mask feature of Figure 2; 

Figure 5 is a calculated starting gray scale iUurnination configuration (J J) for the 
mask feature of Figure 2; 

Figure 6 is a binary representation of the iUurnination configuration of Figure 5; 

Figure 7 shows analysis of a print of the mask feature of Figure 2 printed with an 
annular illuminarion configuration; 

Figure 8 shows analysis of a print of the mask feature of Figure 2 printed with an 
optimized elliptical illumination configuration; 

Figure 9 is a map of the calculated optimized four dimensional illumination 
configuration for the mask feature of Figure 2 scaled to 1 lOnm design rules; 

Figure 1 0 is a calculated starting gray scale illumination configuration for the 
mask feature of Figure 2 scaled to 1 lOnm design rules; 

Figures 1 1 a and 1 lb are binary representations of the illumination configuration 
of Figure 10 with differing values of o; 

Figure 1 2 is an example of a mask pattern with critical gates and cells indicated; 

Figure 13 is the mask pattern of Figure 12 with assist features added to reduce 
the number of pitches in the pattern; 

Figure 14 compares probability density functions of space widths of the mask 
patterns of Figures 12 and 13; and 

Figure 15 is a schematic representation of an apparatus for 
microphotolithography. 

In the various Figures, like parts are identified by like references. 
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1 Abstract 

A method and apparatus for microlithography. The method and apparatus 
include optimizing illumination modes based on characteristics of a specific mask 
pattern. The illumination is optimized by deterrnining an appropriate illumination 
mode based on diffraction orders of the reticle, and the autocorrelation of the 
projection optic. By elimination of parts of the illumination pattern which have no 
influence on modulation, excess DC light can be reduced, thereby improving depth of 
focus. Optimization of mask patterns includes addition of sub-resolution features to 
reduce pitches and discretize the probability density function of the space width. 

2 Representative Drawing ! 

Fig. 1 
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[0 02 6] *7 , *>X©^S;{brit. a#©JHl*4& 
©&3i{fcrtt. 4 yc5cfflEjai»»» ( T C C ) 4*136 ' 

u Tcc©a»7-»;i«gs*sa«toajaitts. tc 

COalfflCC-XsTB. WU*. BorniiiOTo 1 f 
©Principles of Optics (^^© 
US) . 06Ji£. 5 2 8ftt>U5 3 2^~iyicm^ti 

ri>*. cnu, ?iffl(c«t«j*®(cfe#*ti4t©<t-r 

[0 02 7] TCC«. fiPJt><fc#-T3tf? I, 

<&©e assure $> 5 „ n i tc, .3o©a^-<r-5R©5iiE? 

S*> liyi&J, (a. tf) 4*U CCTafc<fct>' 

i8ttiRi98JKttoffi«r**. «©st*©fcaf>K:. J,© 

^S±©^@a r ^S^-raciAS-C#5. * 
3itt©Sf254?T ;teJ:tfM{c**fc a r ©fij.£4 
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[0 02 8] W©RB. ( — m A./P.N A. -nA 
/P,NA) 4tf^<t-r-5SKi><t*K (a. /3 ) 4^ 
-f. A/NA©HSSt-CiE*i{b§nri,^© 
r. K©^Sttl. Or&S,, £ffliJ©Rfi. iBjfiHc&IJ 
Ot*4^-r*S. CtlB. (pA/P.NA. qA/P, 
NA) 44M>£ U-Cl>£„ cn6©St®2offlST 
«. m. n. p. is£V<ittm®<D\mifr*-&-lctm 
U TCCAS±^©<l:^^4^5u (4-D) ©5£-C&S 
io Ck&BfllotPlcteZ. x^[6]©Ii]Sf^--^-tJm*5j:0* 
pfcJc-orSlSft. y^OIilKT^-^-anteic/q 
lcJ:iXm.ZtiZ>. C©StHj©?c&icx*sJ:^y)^S4 

[0 0 2 9 ] 4 -D©j3U<I©.£ (m. n. p. q ) tto 
^T©TCCB, 3-3©^-C©R*5fi^LrCi-S^U4 

#w tea«(oai»T & mmimm&)-e$> a £ see s 

20 D. TCC«8ltStfi«]r*S. a^Kj^^-*->Bi^-c 

3£U> 



[0 03 0] TCC«, iJr*-jf-«I«tt(DOC * tt. TCCtC/f*->©7- 
C) 4*I5£-f SCir, v<$->©^4^tfJ:5{Cffi iKJ:-,Ti?.h5. 
ITSCtmS. S2JCDOCC4*i^-rS„ Cft*30 

OOCC(m. n, p, q) = T(m, n)T ' {- p,-q)TCC{m. n, p, q) 



[0 03 1 ] MK. 3 i>">CC*$W j£,3 tc 



DOCC {m.n.p.q) 



5£<3) 



[ o o 3 2 ] fsjEje^jfiw. smmiamy -t^ztu 

MliUt^ TCCttSlSIS5tBB8frCft4. 
s2*&**©x tepefy©*- -tti-e*!. 5£4 

&foJ:C«E©;^©*-*-^iJB-C*S. 0fl 
mttft©m..,#>?>iE©m..,$-C{c;R.* (-m 



... ^m^+m„.,) . SfcJ;C>'iE©?^©^--*'-^ 
S-C*€>©-C. TCCC^^tt. 2m„.„+ 1 
40 2n.„+ l8r»*2 p..„ + d..,+ lt* 

s©-c. £r©><c*->ig*T*-y-*fm-*-SsB«M: 

m.«. ^mS + m..,©^-, y7?^]tctet»rti 



(7) 



11 



/«« = ««, = /'..or = flOO, 



{ 



r m 2 0 0 2 

12 



3 3 4 8 3 6 



fPM4(l + CT a )"] 

/,«« = "„ax = 9™ = fl°°>\ £ ( 5 ) 

[0 03 3] ^lte<fct>-2£5£3«:ttA-r£<!:. lO*^Kgt»l^^A-SBSHjc>i^©gP^^^SC i 

m a ;> 4 



5£(7) 



•/,(a.P)r(m.«)r(-p-9K 



a + - 



[0 03 4] ato£V0iZVmo£fr<Dm&*2k-?<D XZfr&mL. S8tttg?hS„ ^8^6. ChKl 

H*J.. ,tt. ffi*3-<D®8r*-#- (m. n. P> q) e©£T©^£ 
tcML-CBSB^t^©i'©g(!» (a, f3 ) #ffl(,>eftTI,>* Ci"C. Si^ 

771,1 _ 77.A 



(a. p, '". n > Pi <r) = j s( a ' py{ m - n ) r ' (- p> 



P r NA P„NA 



(m. n. p. q. a. (3) £^fH~& 
pA qX I 



P.N A 



■ P~ 



PyNA 



/(x. 3,) = jf^3z z z z /^lUSH^ (a , p.,„, „, p , q) 

i£(9) 



[0 03 5] B^^^tc^^-c*^^^. J.. t «6yC7C 

[0 03 6] m+p*5j:f>'n + qk:o^-C^}$?:B(5 
CiT. S>f>fltaSl (x. y) m+p = n 50 



-C. m+p=n + q = 0i55 Bf HJft?t*;cDSP#S:|&£T 



13 

mm*- if-) K^fe-r*. c©4swgw»*j. .,.,.* 

rj = m + p => p = tj-m 
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* <fci>¥.&,, S;i ofc-Jren l £ i (x, y) &c-x>-c©5$; 

0) 



f = «+?=> 9 = 



i(n) 



.:/,_ «V,_. i&±<\ 4tA 



y w ,(a.(J.«.n.n-i».6-») ltd 2) 



[0 03 7]Sl2m J.,,.,.11, SlOteiyi ^Aticir, J..,.,.«Al 3(c«t»«J:9«:Ct-rc 
lfC«£o - r^ifc£^&t,;fc&. mijJrO'nfC-^-C© J S&g 1 (x, y)&. 1 4 CCteW-SJ: 5 CC 



5£( 1 4) 



[0 0 3 8 ] MfU„«.,.tt. £0$?* , 

.„.-,„«, ShJAt^-^-T (m. n) (Ci-jtl&ff 
WStl4©t?. ^c#f @5f^-^-«. Stf»fCS{J-r 30 

[0 0 3 9 ] #S©^*->tt^T^ttJ$H>i;W© 

0. J.„., B (a, 13, 77 = 0. 1 = 0) ?;ig#t5C 
iT$ii!)6nS„ SI 5W 77 = OMyi = 0©t§ 

^*istt^gi3nt:*j6-r. j.. t . J0 (a. 0, 



I 



S ~ — — max 



[0 0 4 1] JHWJBttfcj:jyO<*->H3li*3*i4© 
-C. ftiS&ffilE (OPC) *^MT*i0Sf*-y-«: 

JBtt J « . , ** «fc tf'-c * - >{c»T 5KM*<SHS*tft>tt 
£©«fltSie ( D O F ) . 5 A >£siSP ( E O 



r=0. £ = 0) i20©^-^-^/c«DC^% 
S»5C§F-^LAj;t,>BSBflF»i©*(C«fcor. DC*© 

[0 04 0 ] CCDtclb. J ...KiSIMlDCJli 
3W«#-i-S. aJ,.,tSKt. JBWHSWD £©»»*! 
A 1 ) ^tC4is?*5. 



5) 



©i^£-5x.£©{;t. OPC7^-7 ; ttW<, 
<tO-C©./-c#->-C£>3©r . J ,.,tt«iBttlPJJWJRB8 

©&®{t©;rc«>tc&fcil< iRjgT*C £«:&.&„ 
[0 04 2]WSMJ,. t tt. 0ttC»bl©*5H 

J:-3^TCiAi-Ct5. H5f7fe^^ (DOE) (CJ: 
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[00431 W1: 3fcK«l,fcJ, 10 
«>©8£ffi£> ft^g^gi^-^SSfflL-fc. 150 
nm©;<£->=& 1 3 Onmtei^l 1 0 nmf8:ff*i£ 

(cffl/M,, (na) *to. s©*^^- r> r 

1 3 0 nmaS!«-*i>©^'l«^'i'->*^-r o 
[00441 I3K, C<D-??.i> ■ 7 ^-?t©@W 

- ssiivr. @3-c«, s^©^--y-ia. 

(o, o) y-*fcttDcw*#-c*a. tsfttcs 

fctF^L/c^-^-li, (±2, 0) *-^--C*»3. 
^*s^<*->«:t}W-58S^r[fi]©n^=&^-r. ffi© 20 
SM&3f--#-B (+ 1, ±l)t?*«3. f'JTKm® 

-?-KJ;^tDOCCK:fct« fi#tf W0&#^fc 

[0 04 5] H3(C*j{*4@Jra(— ^-»»T (m. 
n) %S13«:fUl/t, J„ t - 1D £#m|-r£C£#T 

CtiZmilCm-f. a4*>6t>*'5J:5«: > J 30 
...-znJcarrsS*©^-^*, (77 = 0. € = 0)*- 
£--C*-5. (0, 0) *S«{cgF#ti-r. 
DOF*a<-Ti. 5*1 5*^-i"J:^K, C©(0, 

o) ^si-flKiiij,. sc£**c 

(0, 0) *-#-£#ItL&t,>£. Mckoym** 

« (77 = 12, € = o) @Sr*-y-r*»). cnu. 

x^lfij{C?Q-3/c^lBl7-<>©^fi£*^-r. frr> 

nm? j >©*BP*ffl5e^-saB©jss»«. < t? = + 1 . 

€ = ± 1 > @»r*-y— C*S. (0, ±2) 08?*- 
#-tt-?>-¥>^ cnj:9Attl,^-?-«, U> 40 
x©77 = o*jj;iy t € = ±2©^c<:*}t,>r*g-^-r5. * 

-3. DOCCOfSB, HM^/££&#^ £*:#«:£'© 

a, 

[0 0 4 6 ] 5U 5£ffll,>T. 1 3 0 n mlSftSiH©*! 

^li^N- $ - > (omm v- 1 a c t #-e# a. 
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j£-ra. cne>©tiR^-r!i<-^s^cftjxr. MWV- 

i*©*,|>«. H^fi£(cA^<«F#r-2). ±i2©<fc*> 
{c. BSB^Oi^ti. yu- • ^-iT-- ;U*fc«2fil«SB§7' 

[0 04 7]ffll>5ggKj£t:t, ^U- • X£--;l/JJS 

B. IWWpIfie&B89H3*af*3Bl|eO. J8iJBJBtR©d>fc < £ 
&J?T-%-©3B#K:-X>-C. 0t5i£> 1 ^©JEJUbl"*;!/* 

«. .^^^-rAtCteWSIalSf^^RT- (DOE) 

B. 05fCiSf Jr^tc^-r-SCi^s-C^-S.. brt'Uft 

#6. m^±^LtB$n05(cm6na^@fx/^4'©— 
*fi«7^ju*ri*!eofc«tc»*sn*. Sot, m 

fe©i-T-2>„ 

[004 8] 2ffiffiR!»7fctt*fl|l>S»^. Jffi 
^S©3SS«:2fil©* (Ogfcttl ) rtSnJ&g-cfcSJg 
ff§BJfl5tt©S*5C 0 £ fcB 1 ©ffi^rfd iJStSfcS) 

B. 0. 8©Hfl*»lR0fci*^. 0. 8*H*..&jrai!IJ 

£tg©5£afii« itc^jij±if en. o. 8*j?s©«iBo 

[0049)^2 : 2ffi©#tt(C^U- • Xtr-jv^m 

i>r. o. 8 8©s^f*Scr*ssb-c. |B]Dn^755 
[0050] yicti-r . 0 6 ©sas-fbs n/cMBj^tt©^ 

eg*. NA = 0. 8*S«tCXA = 2 4 8 nm<D^.f -jV • 

r> h* ■ + > • «; ^^g±©2ffl'7^i'tc 
ot,»r-> JjIz-k,, StfcffS0J©~> s ^ u- h u/ctt 

o. 7*®^7*c©r-. -<t> hJb (p^i) j*^uy^ h • 
■t^^frffli-^c. c©*-r;U"Ctt, ui?xKi. ®Jf^ 
n=l. 7 6-jO. 0 116) ^Wt5$-/7 - O4 0 
Onm©il$-C$>0. n=l. 5 7 7 -j 3. 5 8 8© 
* l Ji"i3>ilo±©n=l. 45- jO. 
Z>ffl<D$-< y'(D6 6 nm©±ic*5. S7*3<fcO*8tc. 
J«tSM?g (a,„ = 0. 5 8*sJ;C>*a out = 0 . 8 8) ft 

6 vicmm t l tcmwmm. ( a „ . « = o . 8 8) ©m* 
*^-n-en^„ H7fcj;j>*8©^tctjt,ir. ^giM 
ilit© «: *j W 5 WM©$S*te J:0'F-jW>(D^5 

jlu— >a>s^?:^-r 0 cn6o@T«, wy,v^ 

Bossung^Pj hB££tmU $Sm<t LT?f 6n 
^iJslw*. ^<B©^Stc-3)^r. go&f VL^^Xif 
*jyVc-$Z>„ commit. DOF?:l$oaxiLtig 
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v Ft?»aofc«t5*a«attBiia (&®&jsi*fg) -e© 

[005 1 1 17KB, m$Sm (cr,. = 0 . 58fcJ; 
tfa.^O. 8 8) *ffll>fc0. 8©NA-C©2fBv io 

• 7 -< -^-v Ceo t,>T. 1 3 0 n mizfiS^OtxA 
*5|g53-gt>'-;3f->©>'5iU-->a>*gS ; &^-r„ C© 
flStRiSJEtt. -0. 4nmH0. 0jimIA*r, *tj 
0. 4um(DDOF*f1-5. U^h©3>F7^b 

g©3>h5^hTH v^d7^Mlf^:7T f£ (ME 
EF) , gim-??^-- K (EL) #i/h;*<,». 

(EOL) ©5a,ffi#;fa2 0 nmT2ibZ>tt%:7nl,X.\,>2> 20 
Chtt, 1 3 0 n mm\W^^C'O^Xt>'fMC : y A 

ft*tt&Hlt/h L- tttf 4 £ . L It 5 A > #fl&© 7 * - 

5=-+ta5^-rssn*i*5©-c. 7-o©5ss{ifett-9> 

[0052]i8TH. 1 3 0 n mRttSm©n^5M 
^yS^->{co(,iT©^5 3. U-i' 3 0. 
8 ©NAT. 0 6©S®{t2fiiBp.^«*ffl^-r. 2ffi 
v^i? • ^ .< — ?i-(co(,ir0^brt,>4„ 30 
J&Ktt. -0. 45/im«+0. 15/imi*JT. 
iftO. 6 mn©DOF£W-f4. 13 8 ©{RES^H 7 

©fe©tifc«-j-si. jKH^bofcjBMjBtttt. attain 

BtKK-^TCME E F#/M < . JgffiftJSIJHJffcK©*: 
^©gtti^^jL- F^tnc <t*m^4„ C©S 
3S{bra?BBt£©SiJ©*J<£tt. S^JPimclt^-r. v-O 

x^^wms.. zvwmiwRteWfi. ^*->_t©-5 40 

-f>*i£g-r4C£fc< EOL*iliJ#nJfigr&4C££ 
Chit «fc9*MftSI+g*g©lli/J>©fc#KjffSP 
£-C$>4„ 

[005 3 ] #13: 0 7fcJ:t>*8K:fcW4 2ffivx^ 
(BIM) {COti-C©^*. ?OAUX-v^^ (C 
LM) tCO(,ir©->5 * ^3 >»»iJt«LAi. S 

^x.ma<D^m-c. y7i-)iT^5ii/->3> 
nSDOF©5Sc#*>e.+^/S.lLm?rSW4± ^(C. 50 
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(0. 0) *-y-©**s2*t'J-'5. ">S^U-'>3 
>*>e>©HK*S^tCj:o-C, (0, 0) +-^-©3t© 

-rft«&c>„ COWrtt, A/[NA (l+a..,) ] 
*iffi©f e^T?. SiS^(C5 -f >(Cf •( if 'J > y£^T 
->fc„ L-#>L-&*S£>. •r-i1f'J>y©7 r *-7-^ -1f-f 
SSSfr D O F *s J: tf-'-f * - >&^g©fc fete . 
(0. 0) *-#-©3fc©fi&®»ft-r-&«fc5K:3!gL, 

cLM©fc&©^5 a > 3 >$s 

•*^^AT155nmTifeofc(77. 5nm<D?V 
A-7-^7>K) . CCf^rB. (0. ±1)*- 
#-#&Ui>£,2f&cA4C £*siatia3»W 6>ti5. L*-L 

feRH->-Jl/«cio-CDOF*SA6C-rs«fc5K:|S3fSU 

[00 54] f(4 : 1 3 0 n mt3:H«©lCi:-oiiT© 
jl U--> 3 155nm«-7|-->'f 

-cmfnhtiL. 0. 8©NAte«fctmtKIB58 (a ln = 0. 
5 8 teckO'a out = 0 . 88) KJ;0. A = 2 4 8 n m 
©SSST?. CLM^aifc. C©iS«^^©CLM 
tt. DOFtMO. 5 umV$>-otc (- 0. 4(imM*> 
6+0. 1 ymj£*> . ^t«RS9§©C LMtt, mViMW 
©B I MtCtb^T. DOF^A#<. -±M^mCXzi 
>h7Xh*5fnTl'fc. ChS, CLM©tt#6*SBI 
Mv^i/«t»3*>^*aTt,^cC£^:7^-r„ h -j I; ■ 
i/SaU-i/H^SIJ, CLM(CJ:6EOLttte*5 
SlfcfltyCttB I MiCj:SEOL^J:0feMh.-CCi4C 
i. te<tc;, C LMfJB I Mtctb^T^Vi? ^ h • rh — 
Jl/©7 >f ^ > ^il?* «fc O+^ic^-T S C £ fcT'i* 

[005 51 W5: 1 3 0 n mtiA,ifimftM'H> -> 
^SSJliCOI,>r©^5 a >3 >^*?r. 0. 8©N 
Ate<fcCXS6tC7n^jS^btiny-l'd<-JU©A = 2 4 8 
nmOgSicrJUtSSLfc. cn6©*S**. 155 
nmJ^—y b — > • f^5^*iJ:i>*5 0%©fa-f< • 
If ^Jb* M-f Z> wi <DVtC fflt, > tc C L M U f- f £ HI— 
(DXsmh&m^-C^* U- h C©^3l^bMHJ 
^fl7EbfcCLM«, 0. 7/im©D0F (-0. 
5/im*6 + 0. 2iLtm)?rWb. 4 0 %©Bft#-C£) ^, 
/c„ B o s s u n g Va v hit. mM£8&&1)W0 . 2 

i-e&-&c<t£inL-ci,>/t 0 jtra^.r. jE5i^i}i©7^# 
S£&s«fc i 5K:U'3 ; 5';i/*gSSSu. BEDIMS 6 
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CiOTS, CD©i4j— tt4U5T5Ci*i-C^5Ci 

-=e-r.>McS-5< OPCB. EOL£iKt*iIE-r£C 
[ 0 0 5 6 ] {SJ6 : llOn mSff^O^gt^iCO 10 

■y- > ^ y > if&mit-r *> fcfefc . j „ o , - , „ *h 9 «c m 

mh. (n=m+p) ^TM^ftK:. y*- 

¥- (£ = n + q) *mm.-ftfaULm-$ „ 13 0nmlS:ti- 
mm<<C^-C(Dm4 i«. ll lCl lOnmglt 

^cc*hrss*©«-^;j\ cn = o. e = o)*-* 

CO(0. 0) *-#-©3fti. DOF«Ci 

Bfc*3*l4. gfc. @9«. (±2. 0) tf-y-Tfc 20 
< (+ 1 . +1) JHEJBtt©*aft«C*fL 
X&k<D&^ttj:Z>Zt*m-?. CftU. NA=0. 8 
© 2 4 8 n mOgltlJ 1 1 0 n ras&H-SiittaiSffnr 
?-5it>5<^Ofc»-C*0. C©fS?&£*iij£-r£/c 
awcte. t>r*»{cm»NA3Wff*o»<». #Hiif§?i!§*8i3£ 
t^fc*iClt^t5t-f'-H. (±2,0)*- 
#-"C*-5„ 0#>bft#£. (±2. 0) 
fiB^JBttOS^gBtc**) (0. 8<ct<1. 0). C 
tiii. cr# CCi&gtefcW-S 1 lOnmg 

9l«aiOS6tt*gc#?Iffi-c*4C4*mly-CI,>4. 30 
[0 0 5 7 ] SCI 5*J<tt;09 ©*£3l£fflOT. 110 
tC. 1 1 Onmn^a53-gUl«C-3l>-CCDSji'fbfiS^ 

Ci*^0"Ct»S. H 1 laic. £©flP«KtK©l^>© 
?Ifie«:iaM8»4Hw-r4. 2 4 8 n mSISfcfliHTK 

fca&JCtt. mi 1 bKm-r^^JC. cr£l. Oil. ^ 
SH-fe** (a©'J>4fltB.*0. 2) 4«-J-4RH93JB« 
£fflt->-5,, 40 
[0 058] * 7- -i tj)\,tH*)\, 
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k *) i&mit-r zm— <Dm.wBVt&£j$& s c t aj-e * 

§1 2(c, S"j -r -r - y-Fg„ g,. 
<fct>*i> 'J • -fe-U-c c?r*tr&i]!&©10!J£jrj 

w-^-^Hiticim, txicm^tcmwiz 

[0059] mw/^z-xDmmtmz&mt-r&ffl 
©#&«. i«a.'<-K:«fcor/'t*->iasf«:*M-r4c 

<bf*-E>. IfcSL'f-W:. AS I C*/c«i&ffii9:it«:oii 

0. 6 H/NAOi'^B 
S„ 01 3r». ^©f5tgL^-*Jjn^.SCi-CLll 2 

^•c. mmww.zm.mt-rzzLifiX'Zh. x^x\ 

HrX • ■J-fVF^ttfitl*. »SL^-*«r *«:<!:*:<* 
mitLtcmW&Vi<D7''a-bx ■ <v 4> K-J t tb$j?-T€>C 

t5tSL^-i^^bL/ctt«-rn.BSBJ (OA I) t 
©ffl^fe-a-lct. ft^©Bl^D0F-7"P-feX • -W>F 

[0 0 6 0] MW^©SiS{b«:iaS-rS/c«>©SiJ©fti 
±fct. ^ifl (SW) ©#SE(C»-^l>fctSt£L^-©lSa 
«cJ:64>©r&So ffcSL^-tt. )U-Mcmn< OPC 

5*^©3lffl«©tt*a?aEIJB» ( p d f > zntiarzc 

-fb-TSCi^-C^S. mm.<Dy -<>*j«fcCf*^P©7 Y> 
A5te|ig-c&Siffi^-r-5i. ll»f*-y-T (m, n) 

O'nCBilifciLTmt^-y-^lf^B-r.E),, CCfwti 

, = SW, + wfc<fcCfP y = SW, + w«. ^fhxteJ: 
O'y^ffiHCfcWS f *X&Z>. 
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P I NA 



PjNA 



P. N/1 



6) 



[0 0 6 1 ] 5$1 7(i. 4ooS;©^f?'i-C^>f3, &^© 10* fc^D-'m* 0 , n^Otibi. 
JB&t . m = n = 0 . m = 0 . n * 0 , m * 0 . n = 0 . * 

... 2 



T (m . n ) = 



[, _ 

<*^(^)*(-^-)(^)-(-7r) 
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